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Abstract. Barchans are crescent-shaped sand dunes forming in aride regions 
with unidirectional wind and limited sand supply. We report analytical and 
numerical results for dune shapes under different environmental conditions as 
obtained from the so-called 'minimal model' of aeolian sand dunes. The profiles 
of longitudinal vertical slices (i.e. along the wind direction) are analyzed as a 
function of wind speed and sand supply. Shape transitions can be induced by 
changes of mass, wind speed and sand supply. Within a minimal extension of the 
model to the transverse direction the scale-invariant profile of transverse vertical 
cuts can be derived analytically. 
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1. Introduction 

The interaction of turbulent wind with sand is known to create a surprisingly rich 
variety of sand structures in many places on Earth and also on other planets. Among 
the most symmetric structures are isolated heaps and barchan dunes. For a schematic 
sketch that introduces some of the terminology see figure 1. By identifying and 
isolating the mechanisms that are essential for dune formation, a mathematical 
minimal model could recently be formulated [15,16], and was shown to reproduce 
generically field observations reporting systematic shape variations [20]. In particular, 
shape transitions from dunes (with slipface) to smooth heaps as a function of sand 
mass were predicted [1,15,16]. The model adapts, combines, and extends earlier 
theoretical developments in turbulence [12, 26] and aeolian sand transport [19, 22] into 
a consistent mathematical description, which allows to identify a weak spontaneous 
symmetry breaking of the turbulent air flow as the origin of the growth instability 
giving rise to structure formation. It moreover pinpoints the mutual competition 
of this symmetry breaking with saturation transients in the sand flux as the basic 
mechanism responsible for shape selection as a function of size and environmental 
parameters such as density of the air, wind speed or sand supply. The characteristic 
length scale provided by the saturation transients is called the saturation length £ s 
[15,16,22]. Due to the ballistic grain motion, it can be much larger than the grain 
diameter d (« 10~ 4 m), which is the only elementary scale in the problem, and can 
thus compete with the turbulent symmetry breaking to cause scale dependence and 
shape transitions of dunes that are of the order of 10 2 m long. 

Below, we address some consequences of (broken) scale invariance by analyzing 
the predictions of the minimal model for the longitudinal shape of a barchan as a 
function of environmental parameters. The next section reviews the major parameter 
dependencies of £ s and deduces some immediate consequences for shape similarities. 
Section 3 reports a systematic numerical study of the influence of dune size, wind 
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Figure 1. a) Sketch of a barchan dune introducing some terminology, 
b) Longitudinal cut along the symmetry plane. The region of flow separation 
is modelled by a smooth continuation of the dune profile [16] and the shear stress 
on the back of the dune is calculated as for a heap with the shape of the envelope 
(dune and separation bubble) with height H e and characteristic length L e . 



speed and sand supply on the longitudinal shape. Finally, in section 4, we extend the 
minimal model to the transverse direction. We close with a brief conclusion. 

2. Similarity hypothesis 

With only one characteristic length £ s in the problem, one expects that shapes of 
dunes of various size and under various environmental conditions can be related to 
each other by similarity transformations. This immediate consequence of the minimal 
model [15, 16] is supported by careful field measurements [20] and, more recently, 
by laboratory experiments with submarine dunes [7,8,10,14]. If this promising 
evidence could be further validated in future work, submarine structures forming in 
unidrectional flows, which were formerly identified as ripples [3,4] on the basis of 
their superficial similarity to aeolian ripples, would actually have to be regarded as 
proper dunes. This could tremendously facilitate the experimental study of dunes 
and ultimately allow us to draw conlcusions about kilometer-sized dunes on Earth 
and possibly on other planets from analyzing centimeter-sized submarine dunes in 
the lab. However, some caution is needed in postulating simple one-to-one relations 
of similar structures under completely different environmental conditions. Arguably 
the minimal model, which has the form of a concise hydrodynamic description of 
sand transport in turbulent shear flow [15], could well be applicable outside the realm 
of validity of the expressions and values derived for its kinetic coefficients, such as 
4, so that similar dune structures may not always be reliable indicators for similar 
elementary transport processes. 

For any fixed given ratio r/rt of the wind shear stress r to the threshold shear 
stress Tt for the mobilization of grains by the wind, the saturation length is predicted 
[15, 22] to scale linearly in Tt divided by the gravitational acceleration g and the 
atmospheric density p a . According to conventional estimates from simple arguments 
for the aerodynamics on the grain scale, r t = 9p'gd [2], i.e. £ s scales linearly in 
the grain diameter d as required by dimensional analysis and also linearly in p' / p&, 
the immersed over the atmospheric density. The dimensionless (generalised) Shields 
parameter 9 [2, 24] accounts for the complicated interactions of the hopping grains 
with the air and the bed. It usually can be taken to be a phenomenological constant 
for aeolian sand transport on Earth. The density ratio p' ' / 'p a ranges from about 2 • 10 5 
for typical conditions on Mars, through 2.2 • 10 3 on Earth and about 50 on Venus, 
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Figure 2. a) Height- length relations of windward longitudinal steady-state 
profiles for various "masses" (1 — 676 m 2 ) and wind speeds (t° /rt = 1.4, 1.6, 
1.8, 2.0, 2.2, 2.4, 2.6) collapse onto a universal curve if heights and lengths are 
normalized to their values H c , L c at the shape transition. Insets: The windward 
profiles of all dunes and heaps are of the functional form (2). b) The transition 
length L c and the minimum heap length Lo 0.7L C scale linearly in (®. 



to 1.7 for submarine flows and can thus dramatically amplify the effect of broken 
scale invariance beyond the grain scale d. Typical grain sizes found in dunes arc 
estimated to be roughly the same on Mars and Earth [6], and similar (or slightly 
smaller) grain sizes have been employed in laboratory experiments with submarine 
dunes [3,4,10]. One can therefore roughly estimate the ratios of the various Tt as 
Mars/Earth/ Venus/ water w 10/1/1/0.5, respectively, if potential differences in the 
Shields parameter are disregarded [13]. Accordingly, we expect £ s ratios (and thus 
size ratios, if size is measured relative to the corresponding grain size) of 

Mars/Earth/ Venus/ Waters 10 5 /10 3 /10Vl (1) 

for dunes of the same shape, with the same r/r t , and with the same degree of 
saturation of the upwind sand flux. This implies that asymptotically large dunes, 
which should have a universal scale- invariant shape, are most easily produced in water 
[3,4,10]. The similarity rule (1) is well supported by a comparison of aeolian and 
submarine dunes [10,14,20]. However, there may be some evidence [5,25] that it 
overestimates the scale factor between dunes on Earth and on Mars. This could be a 
consequence of the considerable differences in atmospheric densities and corresponding 
differences in the dominant sand transport mechanisms and Shields parameters on 
both planets [18]. 



3. Longitudinal profiles 

To make the content of the similarity hypothesis more palpable, we report here results 
for the longitudinal profiles under various wind and influx conditions obtained by 
numerical solution of the minimal model. To simplify the notation we introduce 
the reference shear stress r° over the flat bed and corresponding notation for the 
saturation length £° = £ s (t°), and the saturated sand flux q® = q s (r ). Figure 2a) 
combines steady-state solutions (influx = outflux) for various sizes and wind speeds. 
The insets demonstrate exemplarily that all windward profiles are well represented by 

h{x) « H c cos" (x/L a ) , L a = L c / arccos(2- 1/Q ) (2) 
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Figure 3. Longitudinal vertical profiles of given mass (225 m 2 ): a) for varying 
wind shear stress (steady-state profiles); b) for given reference shear stress 
t° /ft = 2.2 but varying influx. Again, all windward profiles obey equation (2). 



with a w 3.0 for dunes and a w 1.8 for heaps (H c , L c are defined in figure lb). Such 
collapse of the windward profiles is strongly supported by field data [1, 20]. Moreover, 
the height-length relations for diverse masses and wind speeds can be superimposed 
on a single master curve by rescaling lengths and heights by their values H Cl L c at 
the shape transition. The data points at the lower left indicate that in units of £® 
all heaps have roughly the same length L < L < L c « 1.5L , which is of the order 
of the minimum unstable wavelength of a flat sand bed [1] . The latter therefore also 
defines a minimum size for (steady-state) isolated heaps migrating over plane bedrock. 
Figure 2b) demonstrates that both L and L c are linearly related to 

The dependence of the critical height H c and critical aspect ratio H c /L c on 
the wind strength can be estimated analytically by observing that the flux over a 
longitudinal slice has to vanish upon creation of a slipfacc. This requires that the 
minimum shear stress at the feet of the heap, which deviates from r° by a factor 
proportional to H/L, equals r t if H = H Cl L = L c . Or, H c /L c oc 1 — Tt/r°, which is 
in good agreement with our numerical data (not shown) and mimicked by the aspect 
ratios of slices of fixed mass (figure 3a). By a similar argument [15] one can explain 
the affine dependence of the aspect ratio of a slice of given sand mass under a given 
wind speed on the influx (figure 3b). Note that contrary to the profiles in figures 2a, 
3a, the profiles in figure 3b are not in the steady state, but are slowly growing in time. 



4. Transverse profiles 

So far, we have discussed the profiles of longitudinal cuts through the dune, exclusively. 
To understand the transverse dune profile, the weak lateral coupling of adjacent 
longitudinal slices obviously has to be taken into account. In the spirit of a minimal 
description^:, we regard as the conceptionally most important contributions lateral 
grain scattering in the hopping process and irregular lateral perturbations due to 
turbulence. Both effects already exist on a flat bed and give rise to lateral diffusion 
[17]. Though the neglected systematic slope effectsj cannot a priori be regarded as 

| In the minimal model [15], sand transport is modelled as on a flat sand bed, the topography enters 
only via its effect on the wind speed. Direct slope effects onto the flux (difference between uphill or 
downhill hopping, lateral deflection) and onto the local wind direction (wind deflection) are neglected. 
Due to the relatively mild slopes on the windward side, this is a good first approximation. 
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small relative to this lateral diffusion, we prefer to neglect them here to keep the 
model 'minimal' and analytically tractable. Some investigations of how they affect 
the numerical solutions of the minimal model can be found in the literature [9, 21, 23]. 

Within the mentioned restrictions, we make the phenomenologically plausible 
assumption that the lateral scattering is roughly proportional to the longitudinal flux 
<7U itself, so that a net lateral flux will only occcur in response to lateral gradients 
in gy, i.e. q± rts —£±dq\\/dy [17]§. The new phenomenological length scale £± can be 
expected to be a complicated function of the splash process, but will be proportional 
to (though supposedly numerically substantially smaller than) the average saltation 
length, i.e. the effective characteristic hopping length of the grains, which is predicted 
to be relatively weakly sensitive to changes in wind speed [15, 22] for the same reason 
as the grain speed. Therefore, £± will be treated as a constant phenomenological 
coefficient for the following considerations. We identify the spatial changes 

E y = dq ± /dy^-£ ± d 2 qil /dy 2 (3) 

in the lateral flux as a lateral erosion. From this we can deduce the transverse profile 
if we can relate the longitudinal flux to the height profile. Bagnold noticed [2] that 
for a steady-state dune with migration speed v, the ratio q^(y)/h (y) of the sand flux 
over the brink to the height of the brink must be the same for any longitudinal slice 
h y (x) with slipface, or rather, q^(y) = p s vh h (y), where p s is the density of the sand 
bed. Given that practically no sand is lost over the slipface, this immediately follows 
from mass conservation. Extending Bagnold's relation to regions somewhat upwind 
of the brink, we make the assumption that on vertical transverse slices h x (y) (x fixed) 
the longitudinal flux can still approximately be parametrized by 

q\\{y) w p s vh x {y) . (4) 

Our conclusions would not change, though, if we allowed for additional terms a x y 2 , 
b x y 4 with some x— dependent coefficients a x , b x on the right hand side of this relation 
and/or terms a x , b x y 2 on the right hand side of equation (3). 

We can now derive the profile of transverse slices h x (y) close to the crest of a 
large steady-state barchan dune from the requirement that the total lateral erosion 
and deposition (hence the square of E y ) be minimized for a given area of that slice, 

S J dy [(d 2 q\\/dy 2 ) 2 + nh x (y)] = , h x (±Y) = . (5) 

Here Y(x) denotes the windward base-profile of the dune. Using equation (4), the 
solution for the transverse height profile is found to be a quartic polynomial that 
degenerates to a parabola 

h x (y) oc 1 - (y/Y) 2 , (6) 

in agreement with observations, if we require the slopes dh x (y)/dy\Y at the edges 
to be finite [11,20].|| This boundary condition guarantees that the otherwise eroded 

§ Because of the weak dependence of the average grain speed on the wind speed as a consequence 
of the non-trivial feedback of the mobilized grains onto the wind [15,22], this ansatz is essentially 
equivalent to a similar approach in reference [23] involving gradients in the density of saltating grains. 
Compared to the latter it represents a considerable simplification, because it does not oblige us to 
introduce separate dynamic equations for the density and velocity of the mobilized grains, respectively. 
|| The remark after equation (4) gives some hint as to why some of the more complicated relations 
between the transverse and longitudinal flux proposed in the literature [9,21,23] do not give rise to 
any significant numerical deviations of the transverse profile from the present result. Relaxing the 
steady-state assumption should also be uncritical. 
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heap-like horns of the dune can be constantly rebuilt by the lateral flux. Since the 
total flux leaving the horns of width w is of the order of wq® , conservation of horn 
mass requires w{qg/2 — q{) = J dx q± h , where q^ h is the x— dependent lateral flux onto 
a horn, i.e., q± taken at y = ±(W — w). For large dunes, w <C W and the integral is 
estimated as 2p s vL£±H/W. Using v oc L^ 1 [14, 15], we conclude that the transverse 
aspect ratio H/W should, for given wind and influx conditions, not depend on dune 
size if the horn width w does not, in accord with field observations [11, 20]. 

5. Conclusions 

We have presented evidence derived from the minimal model of aeolian sand dunes 
for similarity rules relating the profiles of dunes of different mass and under different 
environmental conditions. Many of the predictions are strongly supported by field 
measurements and by recent laboratory experiments. Some caution is needed, 
however, when applying the model to extreme conditions, e.g. to dunes on Mars. The 
general structure of the (hydrodynamic-type) model equations is expected to pertain 
to situations where particular predictions for their kinetic coefficients cease to apply. 
Our analysis suggests that both longitudinal and transverse vertical cuts through a 
barchan dune have universal profiles independent of dune mass, wind speed, and sand 
supply, because the scale invariance imposed by the turbulent wind is only weakly 
broken by the saturation transients in the sand flux. In general, the aspect ratios of 
the profiles change with these parameters, however. For the longitudinal cuts, the 
dependencies could be rationalized very well on the basis of the original dynamic 
equations of the minimal model, while a minimal, analytically tractable extension of 
the model was proposed for the transverse direction. 
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